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Abstract A transport theory for atomic matter waves 
in low-dimensional waveguides is outlined. The ther- 
mal fluctuation spectrum of magnetic near fields leaking 
out of metallic microstructures is estimated. The corre- 
sponding scattering rate for paramagnetic atoms turns 
out to be quite large in micrometer-sized waveguides (ap- 
prox. 100/s). Analytical estimates for the heating and 
decoherence of a cold atom cloud are given. We finally 
discuss numerical and analytical results for the scatter- 
ing from static potential imperfections and the ensuing 
spatial diffusion process. 

PACS; 03.75.-b Matter waves - 32.80.Lg Mechanical effects 
of light on atoms and ions - 05.60.C Classical transport - 
05.40.-a Fluctuation phenomena and noise 



Introduction 

Atom optics, the coherent manipulation of atomic mat- 
ter waves, is currently heading towards miniaturization 
and integration. Atom guides in both one and two di- 
mensions have been demonstrated Recent 
experiments have studied atoms trapped in electromag- 
netic solid-state hybrid surface guides (or "atom chips" ) 
l^'S-U'El- "^^^ near future may be expected to see the 
trapping and manipulation of Bose-condensed atoms in 
such microtraps |lO| . Decoherence is an intriguing issue 
in this context because of the close proximity of the cold 
atom cloud to the macroscopic substrate, typically be- 
ing held at room temperature. Thermal electromagnetic 
fields thus perturb the atoms, leading to heating, trap 
loss and scattering The thermal noise spectra 

are much larger than the blackbody spectrum because 
the characteristic distances are much smaller than the 
photon wavelengths for the relevant transition frequen- 
cies: the atoms are subject to thermal near fields It 
is therefore of much interest to estimate the length and 
time scales over which the transport of matter waves in 
atom chips remains coherent. 



In this paper, we outline a transport theory for mat- 
ter waves in low-dimensional waveguides close to metallic 
microstructures. Such structures are used in current ex- 
periments to generate magnetic fields or to refiect optical 
fields. We focus on paramagnetic atoms and their per- 
turbation by the thermally excited magnetic near field. 
Scattering rates and noise spectra are computed for a few 
generic geometries: a metallic half-space, a layer, and a 
cylindrical wire. A self-consistent transport theory for 
the atomic Wigner distribution is formulated for a di- 
lute, non-condensed cloud. In the case of a white noise 
spectrum (that turns out to be an excellent approxima- 
tion for the magnetic near field), the transport equation 
is solved analytically, and the atomic decoherence rate is 
determined. We finally present analytical and numerical 
calculations for a waveguide with static roughness (due 
to imperfections in the guiding field, for example). 

1 Model 

1.1 Perturbation of a trapped spin 

We start from an atom trapped in a wave guide poten- 
tial that restricts its motion to one or two dimensions. 
We suppose that the transverse degrees of freedom are 
'frozen out', the atom being in the transverse ground 
state. In the remaining directions, we assume a free mo- 
tion. Examples of such waveguides are linear magnetic 
quadrupole guides formed by currents in one or several 
parallel wires or planar guides mounted above 

a surface 'coated' with repulsive optical ||] or magnetic 
js) fields. Trapped paramagnetic atoms couple to ther- 
mal fluctuations of the magnetic field and will scatter 
from these at a rate given by Fermi's Golden Rule 
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where Bij^uJa) is the magnetic cross correlation tensor 
(spectral density) and huin is the energy difference be- 
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tween the initial and final states |i), |f). The correlation 
tensor Bij{uj) is calculated in subsection |l.2| , where we 
shall see that it depends on the distance to the metallic 
microstructures. Its frequency dependence is quite weak, 
however, because the relevant transition frequencies are 
at most in the GHz range. Magnetic noise at these fre- 
quencies drives spin- flip transitions that, in a magneto- 
static trap, put the atom onto a non-trapping potential 
surface. The concomitant loss rate has been discussed 
previously . We focus here on the scattering between 
quasi-plane waves in the same spin state that remain 
trapped in the waveguide. These transitions correspond 
to much lower frequencies of the order of the tempera- 
ture of the atomic cloud (typically in the kHz range). 
The corresponding photon wavelengths are much larger 
than the characteristic size of the waveguide structures. 
For this reason, we may adopt the magnetostatic ap- 
proximation in order to compute the thermal magnetic 
field spectral density. 



(We have taken the specific resistance of copper gcu = 
1.7xlO"SrJcm[|l).) The geometric tensor Y,j (|^ has 
dimension (1 /length), and its magnitude is the inverse 
of the characteristic scale of the waveguide. As a conse- 
quence, the scattering rate (^ is quite large for a typical 
micrometer-size waveguide. 

1.2.S Typical geometries. We now give the results of 
the evaluation of the geometric tensor ( ^0| ) for three 
generic geometries. For a metallic half-space, we find 
that 



where the diagonal tensor tij has elements (|, Ij,)- Note 
the long range 1/z dependence on distance z 

For a metallic layer of thickness d above a noncon- 
ducting substrate, we find 



1.2 Magnetic near field 

1.2.1 The method. Fluctuation electrodynamics (or 
source theory) is a convenient tool to compute the spec- 
tral correlation tensor of the electromagnetic field in an 
inhomogeneous, dissipative environment. Basically, one 
introduces a fluctuating current density whose spectral 
density is determined from the imaginary part of the di- 
electric function. The field radiated by these currents is 
computed from the Green tensor for the chosen geom- 
etry It has been shown that this scheme also 
provides a consistent way to quantize the electromag- 
netic field in dispersive and absorbing dielectrics |l6|,|l7|, 
|l|,|l9| (see [pOl for further reference). 



1.2.2 Low-frequency noise. For our present purposes, a 
simplified Green tensor is sufficient: we focus on nonmag- 
netic media and low frequencies where the magnetostatic 
approximation is valid. As outlined in appendix we 
find the following expression for the scattering rate (pi) 



7 = CoYn 



where the prefactor is 
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Here, /x„ is the projection of the magnetic moment along 
a static bias field, \s) is the trapped spin state, and T, 
Q are temperature and specific resistance of the metallic 
microstructures. Finally, y„„ is the matrix element along 
the bias field of a geometric tensor Yij defined in (|30|). 
The scattering rate (||) is of the order of 
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This expression reproduces the result (||) for the metallic 
half space when the distance z from the upper interface 
is much smaller than the layer thickness d. On the other 
hand, at large distances a faster decrease c?/z^ takes over. 

Finally, for a cylindrical wire of radius a, we found 
a cumbersome expression involving an elliptic integral. 
Nevertheless, simple results are obtained for the trace 
of the geometry tensor in the following limits: (1) If the 
distance R from the wire axis is large compared to wire 
radius, 
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Note the even stronger power law decrease jB? to 
leading order compared to the planar geometries. The 
three-term expansion (|^) is reasonably accurate down to 
R > 1.6 a. (2) In the short-distance limit R — a <^ a, we 
recover the flat half-space result (||) 
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In Figure 0, we plot the corresponding scattering rates 
as a function of distance (i.e., z for the planar guides, 
R ~ a for the wire guide). The static bias field at the 
waveguide center is taken perpendicular to the metallic 
surface. The scattering rate for the wire is actually over- 
estimated because the trace of the magnetic correlation 
tensor is used. 
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1.3 Transport equation 

The evolution of a trapped matter wave in the waveg- 
uide is conveniently described in terms of a transport 
equation. This equation aUows to characterize the evo- 
lution of the single-particle spatial density matrix (or 
coherence function) 

p(x;s) = (r(x+is)^(x-is)), (9) 

where the average (. . .) is taken over the spatial and 
temporal fluctuations of a perturbing potential. Instead 
of working with the coherence function (^ , we formulate 
the transport equation for the Wigner transform of the 
density matrix 

w^(x, P) = 1 (^^e^^'^^'X^; (10) 

that may be interpreted as a quasi-probability distribu- 
tion in phase space. Here and in the following, the co- 
ordinates X and p describe the motion in the D — \,2 
dimensional waveguide. 

The quasi-free wave function for the trapped spin 
state \s) in the waveguide is perturbed by the magnetic 
potential 

y(x,t) = (s|Ain|s)B„(x,t), (11) 

where -B„ is again the field component along the static 
bias field. We now have to take into account both energy 
and momentum changes in the scattering process. This 
may be done in terms of a master (or transport equation) 
for the Wigner distribution. It is derived using second 
order perturbation theory with respect to y, assuming 
that V has gaussian statistics and doing a multiple-scale 
expansion of the Bethe-Salpeter equation for the coher- 
ence function |Q. We find: 

(9* + -p- V,, + Fe,f Vp)m/(x,p) = (12) 

j d^p' Sv{p - p'; - Ep,) [iy(x, p') - M^(x, p)] , 

with the de Broglie dispersion relation Ep — p'^ /2m and 
where ^(q; AE) is the spectral density of the perturba- 
tion 

Sviq;AE) = ^J f^{Vi^ + s,t + T)V{^,t)) x 

X e-*(q-^-^^^)/'*. (13) 

The average (. . .) is again taken with respect to the mag- 
netic noise field. We assume for simplicity that the noise 
is statistically stationary in time and along the waveg- 
uide directions. This is actually the case for linear or 
planar waveguides parallel to planar structures and for 
a linear guide parallel to a wire. 

The left hand side of the transport equation ([l^) de- 
scribes the ballistic motion of the atom subject to the 
external (deterministic) force Fext- The right hand side 



describes the scattering off the magnetic field fluctua- 
tions. As a function of the momentum transfer q, e.g., 
the spectral density 5y(q;Zi£') is proportional to the 
spatial Fourier transform of the potential, as to be ex- 
pected from the Born approximation for the scattering 
process p' ^ p = p' 4- q. The transport equation thus 
combines in a self-consistent way ballistic motion and 
scattering processes. 

The frequency dependence of the magnetic field noise 
has already been treated in the previous subsection. We 
saw that the spectral density is flat in the magnctostatic 
approximation. The magnetic perturbation can therefore 
be treated as a white noise. To get the wave vector de- 
pendence, we have to evaluate the spatial Fourier trans- 
form of the two-point correlation function i?ij(xi, X2). 
In the magnctostatic limit, it turns out that above a 
planar metallic layer, the correlation function is simi- 
lar to a lorentzian, with a correlation length Ic parallel 
to the layer given by the distance z perpendicular to 
the layer (see figure ||). More details are discussed else- 
where [|l3|. As a consequence of the finite correlation 
length, the scattering kernel 5(q; AE) in the transport 
equation (^2|) exhibits a cutoff for momentum transfers 
|q| > h/lc ~ h/z. The spatial smoothness of the mag- 
netic field thus limits the possible scattering processes, 
although arbitrary energy transfers are available from 
the white noise perturbation. 

In the following, we assume a constant force Fext 
and focus on two limiting cases of the transport equa- 
tion (|l^): (1) 'Inelastic transport' under the influence of 
white noise perturbations. This is the appropriate limit 
to estimate the detrimental effects of thermal magnetic 
fields. (2) 'Elastic transport' where a time-independent 
perturbation is assumed. This allows to describe in the 
same framework the influence of static roughness in the 
waveguide potentials from which the atoms scatter. In 
optical potentials, roughness of this type occurs when 
the optical fields are scattered from small-scale inhomo- 
geneities in the microstructure p3t . In planar gravito- 
optical traps based on a planar evanescent wave mirror, 
the rough optical potential couples efficiently the atomic 
motion normal and parallel to the mirror . The trans- 
port equation is valid as long as one is interested in the 
evolution on length scales large compared to the correla- 
tion length of the roughness. In rough evanescent fields, 
we find again that this length is typically of order z 
because high spatial frequencies K ^ 2tt/X give rise to 
exponentially damped fields oc exp(— iiTz). 

2 Results 

2.1 Inelastic transport 

2.1.1 Analytic solution of the transport equation. For a 
broad band spectrum of the perturbation, we may ne- 
glect the dependence of 5y(q, Zi£') on AE. The trans- 
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port equation then simphfies because the integra- 
tion over the initial momentum p' is not restricted by 
energy conservation. Taking the Fourier transform of the 
Wigner function W{x,p) with respect to both variables 
X and p (with conjugate variables k and s), it is simple 
to derive the following solution 



W(k, s; t) = Wo(k, s - hkt/m) e-'^oxt-st/?! 



X exp 



-7 / (1 - C(s - hkt'/m))dt' 
Jo 



(14) 



Here, M^o(k, s) is the double Fourier transform of the 
Wigner function at initial time t = 0, and 7 and the 
normalized spatial correlation function C(s) are related 
to the correlation function of the perturbation by 

7C(s) = 1 r^dr {V{s, t) F(0, 0)), C(0) = 1. 
ri J -00 

(15) 

We also note that 7 is the transition rate 7i^t for scat- 
tering processes with zero energy transfer (ton — 0). For 
a wave guide above a metallic surface, the estimate (H) 
is hence applicable for the rate 7. 

From the analytic solution (p^, it is easily checked 
that in the absence of the perturbation, the spatial width 
{6x'^{t)) of a cloud increases balHstically according to 
{Sx'^lt)) = {Sp^{0))t^/m^ where (V(0)) is the initial 
width of the cloud in momentum space (this latter width 
remains constant in this case, of course). 

2.1.2 Spatial decoherence. More interesting informa- 
tion may be obtained for a nonzero scattering rate 7. 
Note that the spatially averaged atomic coherence func- 
tion is given by 



d'°a;p(x;s) = W^(k = 0,s). 



(16) 



The solution (JlJ) therefore implies that the spatial co- 
herence decays exponentially with time: 

r{s; t) = ro{s) exp [-7t(l - C(s)) - iF • s t/h] . (17) 

The decoherence rate depends on the spatial separation 
between the points where the atomic wave function is 
probed, and is given by 7(3) — 7(1 — C(s)). It hence 
saturates to the value 7 at large separations and de- 
creases to zero for s — > 0. The decay of the coherence 
function (17) is illustrated in figure^. One observes that 



at time scales t > I/7, the spatial coherence is reduced 
to a coherence length ^coh ~ Ic- After a few collisions 
with the fluctuating magnetic field, the long-scale coher- 
ence of the atomic wave function is thus lost and persists 
only over scales smaller than the field's correlation length 
(where different points of the wave function 'see' essen- 
tially the same fluctuations). For larger times t ^ I/7, 
decoherence proceeds at a smaller rate that is related to 
momentum diffusion, as we shall see now. 



2. 1.3 Momentum diffusion at long times. The be- 
haviour of the atomic momentum distribution at long 
times may be extracted from an expansion of the an- 



alytic solution (17) for small values of s. Assuming a 
quadratic dependence of the field's correlation function, 
C(s) Ri 1 — s'^/ll, as one would expect for Lorentzian cor- 
relations, we find that the atomic momentum distribu- 
tion is gaussian at long times; it is centered at po-|-Foxt^ 
due to the external force, and its width increases accord- 
ing to a diffusion process in momentum space 



(Sp'it)) « {Sp'm 
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This was to be expected: the atoms perform a random 
walk in momentum space, exchanging a momentum of 
order Ti/lc per scattering time I/7. The momentum dif- 
fusion coefficient Dp = Ti'^/ll that may be read off 
from (|l8|) is consistent with this intuitive interpretation. 
Physically speaking, the atomic cloud is 'heated up' due 
to the scattering from the fluctuating potential. We note 
that the rate of change of the atomic kinetic energy in the 
wave guide plane is the same as the one for the tightly 
bound motion perpendicular to the metallic surface (see 
12 for a calculation of this rate). 

Translating the width of the momentum distribution 
into a spatial coherence length, we find a power-law de- 
cay at long times, ^coh = Icl \fjt- Finally, a similar cal- 
culation yields the width of the atomic cloud in posi- 
tion space: it increases 'super-ballistically' at long times, 
{5x^{t)) (X t^, as a consequence of heating. 



2.2 Elastic transport 

A waveguide potential with static roughness leads to 
energy-conserving scattering. The integration over the 
initial momenta p' in the transport equation (|l^ ) is then 
restricted which complicates the analytic solution. We 
discuss in the following numerical results and an ana- 
lytic solution for vanishing external force. 

2.2.1 Numerical results. The transport equation may 
be solved using a split-step propagation method. We al- 
ternate between ballistic motion and the scattering. In 
one dimension, e.g., the transport equation reads 



P 



{dt + -d,+F,^tdp]W{x,p,t) 
7(p) {W{x,-p,t)-W{x,p,t)), 



(19) 



where the scattering rate 7(p) depends on momentum 
via the wavevector transfer \2p/h\. The solution for pure 
ballistic transport is described by 

W{x,p,t + At)^ (20) 
W{x - ip/m)At + {F,^t/2m)At^,p- F,^tAt,t), 
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while the scattering kernel gives 

/ W{x,p,t + At) 
\W{x,^p,t + At) 

--yip) At ( cosh^At sinh7z4A 
I sinh ^At cosh ^At J 



(21) 



W{x,p,t) 
W{x, ~p, t) 



For sufficiently small time steps At, the combination of 
these two steps gives an efficient way to compute the 
time evaluation of the Wigner distribution. 

We present results for the two cases of an atomic 
sample moving in a wave guide with a static rough po- 
tential with and without an external force applied that 
accelerates the sample. In the first case we observe that 
for short times the average motion is ballistic, which 
is shown in figure IJ: the cloud is accelerated accord- 
ing to {x(t)) = (x(0)) + {p{0)/m)t + (Fcxt/2m)<2. How- 
ever, at longer times a different behaviour takes over 
and the average position shows a linear drift. We also 
observe a monotonic increase with time of the spatial 
width {6x'^{t)) of the atomic cloud. Figure ^ shows that 
this increase is approximately linear, as would be the 
case for a spatial diffusion process. The transient ballis- 
tic regime may be understood qualitatively by inspec- 
tion of the average momentum {p{t)) and momentum 
width {Sp'^(t)), also shown in figure ^ as well as of the 
snapshots of the temporal evolution in phase space (fig- 
ure H): the elastic scattering transfers atoms from neg- 
ative velocities to positive velocities until the average 
velocity vanishes. This is accompanied by an increase 
of the momentum spread. At the end of the transient 
regime, both the momentum spread increase slows down 
and the motion in position space is no longer ballistic. 
Note that the deceleration by the scattering potential 
is faster than expected from the short time acceleration 

(pit)) = (p(0)) + Fexti. 

For the case of a vanishing external force we now 
show that the characteristic behaviour of the atomic 
cloud can also be obtained analytically. The discussion 
is restricted to D = 1 for simplicity. 



2.2.2 Analytical estimate of the diffusion coefficient. In 
the absence of an external force, the momentum p en- 
ters the transport equation only as a parameter, and we 
are left with a coupled system for the Wigner functions 
W{x, ztp, t) [cf. also (|2|)]. Using a spatial Fourier trans- 
form and a temporal Laplace transform (with conjugate 
variables k and C), it is straightforward to obtain the 
following solution 



W{k,±p; = 

(lip) + C - ikp/m)Wo{k, +p) 



(22) 



■-f{p)Wo{k,-p) 



+ 2"f{p)C + k'^p^/nfi 

where WQ{k,+p) is the initial (Fourier transformed) 
Wigner function. The smallest of the denominator's 
roots Ci,2 determines the long-time behaviour of the 



Wigner function. For small values of k (corresponding 
to large spatial scales), it is given by 



k'^p^ 



2to^7(p) 



(23) 



For a Wigner function with gaussian initial data, we thus 
find the following long-time asymptotics: 



t 



{x{t)) = (x(0)) + 
{5x^{t)) = {Sx^G)) 



P 



2m^{p) ' 
pH 



2rn^7(p) 



(24) 
(25) 



The average position thus tends to a constant, that is dis- 
placed by the free flight distance during the scattering 
time 1/7(73). For longer times, the 'forward' and 'back- 
ward' moving Wigner distributions W{x, ±p, t) have the 
same weight, and there is no spatial displacement on 
average. The spatial width increases linearly with time, 
in accordance with a random walk: for each scattering 
time l/7(p), a position step ^ p/'m'y{p) is performed in 
a random direction. 

In figure |^, we compare these results to a numerical 
solution of the transport equation without external force. 
The predictions ( |2^ ) for the average position and ( |25| ) for 
the spatial width show a quite good agreement especially 
for the average position. The transient regime that is re- 
sponsible for equilibration of the velocity classes however 
is not captured by the analytical solution and thus the 
spatial width is only a good estimate for the numerical 
solution, but definitely both exhibit the predicted linear 
behaviour. Figure ^ shows the dynamical evolution of 
the Wigner distribution, illustrating the growth of mo- 
mentum classes with opposite signs. 

The structure of the Wigner distribution in phase 
space when the second momentum class builds up also 
modifies the spatially averaged coherence function (|l^ ) 
as can be seen in figure ^j. Since the scattering conserves 
energy, the momentum width of the cloud and hence 
the spatial coherence length remain constant in time, in 
distinction to the inelastic scattering shown in fig.^ The 
oscillations of the coherence function are a reminiscence 
of the standing wave formed by the two velocity classes. 
Note, however, that this interference is 'classical' in the 
sense that the Wigner function is positive all over phase 
space in this case. 



3 Conclusion 

In this paper, a transport theory for matter wave trans- 
port in low-dimensional waveguides has been formulated 
and solved in a few interesting limiting cases. We have 
shown that the scattering from thermal magnetic near 
fields that 'leak out' of metallic microstructures at room 
temperature limits the coherence of a cold atomic cloud. 
After a few scattering times, the spatial coherence of 
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the ensemble gets reduced to the correlation length of 
the magnetic noise. This length scale is typically com- 
parable to the distance of the waveguide to the metallic 
structures. Our results indicate that decoherence may 
be reduced by working with smaller metallic structures, 
reducing their temperature and their specific conductiv- 
ity. Hopefully, a reasonable compromise between these 
conflicting requirements may be found. 

The decoherence rates obtained in the present frame- 
work are quite large and should give rise to measurable 
effects when interference experiments with trapped mat- 
ter waves are performed. From the theoretical point of 
view, it is imperative to generalize the transport the- 
ory to Bose-condensed samples: the transport equation 
then becomes nonlinear, and the potential fluctuations 
may expel particles out of the condensate. This issue is 
currently under study and results will be reported else- 
where. Finally, it is well known that in one dimension, 
static noise suppresses spatial diffusion, leading to An- 
derson localisation j25j. This effect is beyond the scope 
of the present (semiclassical) theory, and its observation 
in atomic waveguides would be a clear manifestation of 
quantum-mechanical interference in multiple scattering. 
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A Low-frequency magnetic noise 



As outlined in subsec. 1.2, we model the thermal ex- 



citations of the metallic structure by currents with a 
correlation function 



0T(xi;w)jj(x2;u;')) ^ 
4:TrheoLj'^fi{uj) 6ij 5{uj - 



y)Im£(xi)J(xi-X2),(26) 



where n{ijj) is the Bose-Einstein occupation number. The 
current correlation is (5-correlated in space provided the 
dielectric function is local. In the magnetostatic approx- 
imation, the vector potential generated by the currents 
is given by (in SI units) 



(27) 



The cross correlation tensor of the vector potential 
is then simply obtained from the thermal average of 
A*(xi;ti;)Aj(x2; w'). Taking the curl with respect to 



both xi and X2, we get the magnetic cross correlation 
tensor: 



(i3*(xi,u;)i?,(x2,c.')) = 

= 2ti5{uj - tj')By (xi,X2; w), 

Sy(xi,X2;w) = S'g '(uj)- -Y,, 



47rcj/c 



^3^, (Xi -X0,(X2 - X')j 



X2 



^'|3 



(28) 
(29) 
(30) 
(31) 



We have normalized the spectral density to Planck's 
blackbody formula 



S 



(66) 



37reoc5(e^'^/fcBT _ i) 



(32) 



The integration in (31 ) runs over the volume V occupied 
by the thermal currents (where the imaginary part of 
e(x', w) is nonzero). Putting xi — X2 = x and using the 
scattering rate (|^), we find from ( pof ) that the scattering 
rate is of the form (|^). 

For homogeneous metallic structures, the dielectric 
function is e(x',a;) = [/(eoguj) where g is the DC con- 
ductivity. In addition, the relevant frequencies are so low 
that the high-tcmpcrature limit of the Planck function 



is applicable. The prcfactor in (29) then approaches the 
constant value 



U! 



0: S^^'\u;) 



3Ime 
Attlo/c 



(33) 



and we find the prefactor in (|^). We have compared the 
results of the magnetostatic approximation to exact cal- 
culations using the retarded Green function for a layered 
medium. Both calculations agree in the short-distance 
limit where the distances are small compared to both the 
wavelength and the skin depth S = y^2eog/ (ujc^). How- 
ever, the magnetostatic approximation overestimates the 
field components parallel to the layer by a factor of three. 
This is related to the fact that the noise currents ( |2^ ) 
are not divergence-free at the surface. Possible improve- 
ments will be discussed elsewhere in more detail. 
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Fig. 1 Scattering rate in a waveguide at distance z from 
three metallic structures: half-space, thin layer, thin wire. 
The dashed line gives the large-distance limit (|7|) for the wire. 
The metal is copper at room temperature. Both the layer and 
the wire have a thickness 1 ^xn. 



exact C(s) 
asymptotics 




2z 4z 6z 
sepaiation s 



vacuum 
suiface 



Fig. 2 Spatial (normalized) correlation function of the ther- 
mal magnetic near field above a metallic surface at frequency 
uj/2ti = 30 MHz (copper at 300 K). The separation s gives 
the distance between two observations points at the same 
height z above the surface. Dots: exact evaluation, solid 
lines: Eisymptotic expansions in the short-distance regime, 
discussed in jl^]. For lower frequencies, the correlation func- 
tion remains essentially unchanged. 




Fig. 3 Illustration of spatial decoherence in an atomic wave 
guide. The spatially averaged coherence function r(s,t) is 
plotted vs. the separation s for a few times t. Space is scaled 
to the field correlation length Ic and time to the scattering 
time 1/7. A Lorentzian correlation function for the pertur- 
bation is assumed. 





Fig. 4 Numerical calculation of the temporal evolution of 
an atomic cloud in a one-dimensional waveguide. The cloud 
is subject to a positive external force and elastic scattering by 
a static rough potential, with a negative initial average veloc- 
ity. The two lower curves in the left graph show the evolution 
of the average position {x{t)) determined numerically (solid 
line) and the ballistic short time behaviour (dotted line). The 
upper left curve displays the simulated variance {5x'^{t)) of 
the width in position space. The two lower curves on the right 
are the numerical values (solid line) for the average momen- 
tum {p(t)) of the cloud and the according short time ballistic 
acceleration (dotted line). The upper curve on the right illus- 
trates the numerically computed width in momentum space 

The scattering rate is taken as 7(p) = 706 '^''^ with Ic = 
0.1 units, 70 = lunit. 






Fig. 5 Snapshots of four stages of the temporal evolution of 
a Wigner distribution accelerated by an external force. The 
underlying dataset is the same as the one used for figure ^ 




time |arb. u.| 

Fig. 6 A free evolving Wigner distribution subject to elastic 
scattering, the distribution having initially a negative average 
velocity. The two lower curves show the evolution of the av- 
erage position {x{t)) determined numerically (solid line) and 
from the analytical solution (^) (dotted line), whereas the 
upper curves displays the numerical average width {Sx^{t)) of 
the distribution (solid line) and its analytical estimate from 
iM) (dotted fine). 
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Fig. 7 Snapshots of four stages of the temporal evolution 
of a Wigner distribution only subject to static roughness and 
vanishing external force. The underlying dataset is the same 
as the one used for figure ^ and ^. 
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Fig. 8 Illustration of spatial decoherence in an atomic wave 
guide, caused by elastic scattering due to static roughness (no 
external force). The spatially averaged coherence function 
r{s,t) is plotted vs. the separation s for a few times t. Space 
is scaled to the field correlation length Ic and time to the 
scattering time I/70 = 1/7(0). The curves were extracted 
from the dataset underlying also figure pi 



